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EFFECTS OF FLUID PROPERTIES AND GRAVITY LEVEL ON BOILING 
IN THE DISCRETE BUBBLE REGION 
by Thomas H. Cochran and John C. Aydelott 
Lewis Research Center 

SUMMARY 

The effects on boiling of changes in fluid properties and gravity level were investi- 
gated for liquids of high viscosity (a sucrose solution) and low surface tension (an 
ethanol-water solution) in relation to water. The zero-gravity data were obtained by 
allowing the experiment package to free fall in a 2. 25-second drop tower, which permit - 

_Cj 

ted the attainment of less than 10 times Earth gravity. Data taken from high-speed 
motion pictures indicate that viscous effects were negligible during the growth of bubbles. 
However, in zero gravity, the viscous effects were important near bubble separation for 
the sucrose solution. Bubbles generated in the ethanol-water solution had similar maxi- 
mum radii and lifetime characteristics in normal and zero gravities. 


INTRODUCTION 

The development of vehicles for space exploration has prompted investigation of the 
effect of gravity level on the transfer of heat to fluids. A particular area of interest is 
the storage of cryogenic liquids in closed containers over long periods of time, as would 
occur on a coasting space flight. Under such conditions, energy input from solar radia- 
tion will cause the tank pressure to increase, and as a result, the liquid may become 
subcooled (ref. 1). Consequently, a heat-transfer phenomenon of importance in reduced 
gravity is the subcooled nucleate boiling process. 

In a recent study (ref. 2), the authors investigated the effect of subcooling and grav- 
ity level on the boiling characteristics and bubble dynamics of water vapor generated on 
a flat horizontal surface. The results indicate the gravity- independent nature of boiling 
at high subcooling conditions and the importance of a newly defined pressure force on 
bubble dynamics. 

The present work is concerned with the effect of liquid properties and gravity level 


on the dynamics of bubbles formed on a flat horizontal surface. The liquids investigated 
were a 60 percent by weight sugar-water solution, sucrose, and a 10 percent by volume 
ethanol-water solution. In relation to water, the sucrose solution had a greater viscos- 
ity, while the ethanol-water solution was characterized by a reduced surface tension. 
High-speed motion pictures were taken at subcoolings of approximately 6°, 15°, and 

21° F (3. 34°, 8. 34°, and 11. 67° C), in both normal gravity and in an environment of less 
-5 

than 10 times normal gravity. The heat flux selected was in the range defined in nor- 
mal gravity as the discrete bubble region (ref. 3). For such heat fluxes in normal grav- 
ity, the bubbles generally grow singularly on the heated surface, and the heat -transfer 
mechanism is dependent on the disturbance in the liquid caused by the growth and depar- 
ture of the bubbles. 

The photographic data were analyzed by measuring the bubble growth characteristics 
and calculating the forces acting on the bubbles during their growth. 

APPARATUS, EXPERIMENTAL PROCEDURE, AND DATA REDUCTION 

Because the apparatus, experimental procedure, and data reduction procedures are 
presented in detail in reference 2, only a brief description of them is given herein. 

The zero-gravity data were obtained in a drop tower by allowing the experiment 
package to free fall unguided for 85 feet (25. 91 m). A gravity level of less than 10”^ 
times Earth gravity (termed zero gravity in this report) resulted from the package falling 
in a protective air-drag shield. Deceleration occurred after 2. 25 seconds when wooden 
spikes on the drag shield embedded in a box of sand. The experiment package contained 
a boiling apparatus, camera and lighting equipment, power supplies, and associated con- 
trols. Within the boiling apparatus, which is shown in figure 1, were the primary heater 
(a Chromel strip 0. 5 in. wide by 0. 005 in. thick by 0. 5 in. long (1. 27 by 0. 0127 by 
1. 27 cm)), a secondary heater which was used to raise the bulk of the liquid to the satu- 
ration temperature, a thermistor to monitor the bulk temperature, and a thermocouple 
attached to the underside of the heater. The 16-millimeter motion -picture camera pro- 
vided a filming rate of approximately 6500 pictures per second. 

Prior to a data run the glassware was cleaned and the primary heater was polished 
and rinsed with ethanol. The test liquids were deaerated by boiling the sucrose solution 
for approximately 1 hour and the ethanol-water solution for approximately 20 minutes. 

The ethanol-water solution was boiled for a shorter time to maintain the alcohol concen- 
tration as near 10 percent by volume as possible. After filling the boiler and raising the 
apparatus to the top of the drop tower, the liquid was heated with the secondary heater to 
approximately the saturation temperature. Power to this heater was removed, and the 
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primary heater was turned on and set at a power level that initiated boiling. Because the 
heat input by the primary heater was not sufficient to maintain saturation conditions, the 
liquid bulk became subcooled. The power level to the heater was then increased to the 
desired level and the package was dropped. Motion pictures were taken during the last 
1.25 seconds of zero-gravity time. Normal-gravity testing was the same as that just de- 
scribed with the addition that the thermocouple was monitored with a millivolt potentiom- 
eter. 

The bubbles recorded on the film were viewed and measured on a motion analyzer 
that magnified the image eight times. A statistical analysis involved the recording of 
bubble lifetime and maximum radius for as many as 15 bubbles at each test condition. 

As a result of this study, an average bubble at each test condition was selected for a 
force analysis. 

The forces calculated (ref. 2) were 
the buoyancy force 


f b = (V - v b )p, JL 

g c 


( 1 ) 


the surface-tension force 


the pressure force 


the drag force 


F sy ~ <T S 7rD b 


sin cp 



( 2 ) 

(3) 


(4) 


and the dynamic force 


dy 


_ Py /dV dy [ y d 2 y\ 


g„ \dt dt 


dt 


+ F dr + F sy - F B - F P 


(5) 


(All symbols are defined in appendix A. ) A slightly different approach in determining the 
viscous effects on the bubble than that taken in reference 2 was applied and is outlined in 
the section Effect of Viscosity on Bubble Dynamics. Force histories (plots of force 
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against time) were prepared for each bubble so that the effect of fluid property changes 
could be studied. 


DETERMINATION OF FLUID PROPERTIES 

A study was made to determine the pertinent physical properties of the test liquids. 
Surface tension as a function of temperature was calculated for the different liquids and 
is presented in figure 2(a). The method of calculation is given in appendix B. Changes 
in viscosity with temperature were described in reference 4 for the sucrose solution and 
in reference 5 for water. Values for the ethanol-water solution (fig. 2(b)) had to be ob- 
tained from measurements made with a rolling-ball-type viscosimeter. The specific 
gravity of water was taken from reference 6 and that of the ethanol-water solution was 
calculated from the values for water and pure ethanol (refs. 6 and 7). The specific grav- 
ity of the sucrose solution was obtained from both the literature (ref. 8) and experimental 
data, as shown in figure 2(c). 


RESULTS AND DISCUSSION 
Boiling Characteristics 

General . - The temperature and heat -flux measurements taken during the normal 
and zero-gravity tests are presented in table I. Approximate saturation temperatures 
for use in computing the subcoolings were obtained by recording the highest bulk temper- 
ature attained after boiling for 5 to 10 minutes at a full boil. As in reference 2, average 
surface temperatures were tabulated because of the fluctuations. 

Sucrose solution . - The sucrose -solution films indicate that only continuous growth 
bubbles were generated at all test conditions. Coalescence was evident in all the films 
and occurred most frequently at low subcooling in zero gravity. However, coalescence 
never became as frequent as it did under comparable conditions with water because the 
population density of the bubbles in the sucrose solution was much smaller. In refer- 
ence 1, two types of coalescence are described. The first type, which occurred at low 
subcooling in normal and zero gravities, is characterized by two or more adjacent bub- 
bles combining on the surface to form an irregular vapor mass. The second type was 
seen only at low subcooling in zero gravity and was caused by the bubbles moving so 
slowly away from the surface after separation that subsequent bubbles forming on the 
heater surface were eventually absorbed into the rising vapor mass. Both of these types 
of coalescence were observed in the sucrose solution along with a third type, which was 
present in both normal and zero gravities. In this third type of coalescence, an initial 
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bubble formed on the surface and departed; just as this bubble separated, a second one 
was generated from the same nucleation site. Since the growth rate of bubbles is a maxi- 
mum early in their lifetime (ref. 2) and is relatively large for the sucrose solution 
(ref. 9), the second bubble rapidly merged with the underside of the initial bubble, and an 
irregular vapor mass was formed. This process is shown, in a series of photographs, 
in figure 3. Therefore, the primary difference between the second and third types of co- 
alescence is that the cause of the former is the slowness of the bubbles caused by the ab- 
sence of gravity as they depart from the vicinity of the heating surface, while the cause 
of the latter is the high initial growth rate of secondary bubbles. 

The results of the sucrose-solution statistical study indicate the same general trends 
with subcooling as for water. As shown in figure 4, the average maximum radii and life- 
times of both liquids became smaller in both gravity levels as subcooling was increased. 
The effect of gravity level for the sucrose solution, however, was quite different from 
that for water. The sucrose bubbles were, on the average, larger in zero gravity than in 
normal gravity, but the considerable increase in size at low subcooling with a reduction 
in gravity evident in the water data did not occur for sucrose. This effect has also been 
noted by Keshock and Siegel (ref. 9). Examination of average bubble lifetime shows a 
considerable increase at zero gravity compared with normal gravity, particularly at low 
subcooling. Therefore, in a sucrose solution, boiling is gravity dependent close to satu- 
ration conditions, as evidenced by the fact that the bubbles remain on the heater surface 
approximately twice as long in zero gravity as in normal gravity. 

Ethanol-water solution . - The ethanol-water films revealed that the majority of the 
bubbles generated under all test conditions were of the continuous growth type. Occa- 
sional coalescence of adjacent bubbles occurred at low subcooling in normal gravity but 
little coalescence was evident at high subcooling in either normal or zero gravities. 

More coalescence was noticeable at low subcooling in zero gravity; however, it never 
became frequent enough to dominate the phenomenon. Also some oscillators (alternately 
growing and collapsing bubbles) were observed in the low-subcooling zero-gravity test. 
Frequency of bubble generation was greater with the ethanol-water solution than with 
either sucrose or water. The bubble frequency became so great at one site it appeared 
that more than one bubble was being generated at the same time. Keshock and Siegel re- 
ported this phenomenon in reference 9 as a dominant characteristic of boiling pure 
ethanol. The necking down and large contact angles reported in reference 2 as charac- 
teristic of water bubbles at separation were not observed for bubbles generated in the 
ethanol-water solution. 

The statistical study Of the ethanol-water solution revealed the usual changes in bub- 
ble average maximum radius and lifetime with subcooling, as shown in figure 4. How- 
ever, comparison of the normal and zero-gravity data indicates little difference for 
either the average maximum radius or average lifetime curves. In fact, the difference 
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between the normal and zero-gravity results can be attributed to statistical scatter. No 
importance should be attached to the apparent reversal in the curves as compared with 
those for the water data. Care should be taken at this point not to assume gravity inde- 
pendence. Although similarity at both normal and zero gravities is certainly a necessary 
condition for gravity independence, it is not sufficient. In addition, the forces acting on 
the bubbles must be investigated to ascertain the importance of gravity. 


Effect of Viscosity on Bubble Dynamics 

In reference 2, the viscous effects on a bubble were accounted for by a drag force, 
by assuming that liquid flowed around the bubble. After bubble profiles were observed 
(see fig. 5), it was concluded that this assumption was invalid during the growth of a bub- 
ble and that the liquid flow was more closely approximated by a source flow. Therefore, 
viscous effects in the early portion of bubble lifetime can be determined from the equa- 
tions governing liquid motion. For a Newtonian, incompressible liquid with constant vis- 
cosity, and when body forces are neglected the pertinent equations are 

Continuity: 





( 6 ) 


Navier-Stokes: 


Pi 



+ 








2_ ^0 
r 2 00 



(7) 

where r is measured radially away from the nucleation site and 0 is the counterclock- 
wise angular displacement of the r vector from the heater surface. 

If it is assumed that v^ = 0, these equations reduce to 


Continuity: 



( 8 ) 
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Navier -Stokes: 


P 



+ 







(9) 


Substituting equation (8) into equation (9) results in further reduction to 



It is not necessary to solve these equations to ascertain the importance of viscosity. 
Instead, it can be obtained by comparing the magnitude of the viscous and liquid inertia 
terms in equation (10). This comparison was made for a sucrose bubble generated at 
5° F (2. 78° C) subcooling in normal gravity. The terms were evaluated for four values 
of 0 from 10° to 90° and for eight time increments during the growth of the bubble. 

A representative cross section of the results, presented in table n, indicates that the 
absolute magnitude of the viscous terms is always much less than that of the liquid in- 
ertia terms. Therefore, the viscous effects in a direction perpendicular to the heater 
surface can be neglected during the growth of a bubble on a surface. 

For the collapse portion of bubble lifetime, it was apparent (fig. 5) that there was 
some evidence of liquid flow around the bubble. Therefore, for this portion of the bubble 
history, the viscous effects were approximated by the drag formula presented in equa- 
tion (4)’. Force histories for the sucrose bubbles at 15° F (8. 34° C) subcooling in normal 
and zero gravities are presented in figure 6. It is evident that the bubbles were gravity 
dependent as indicated by (1) the large buoyant force in normal gravity and (2) the dissimi- 
larity in the force histories between normal and zero gravity conditions. The normal 
gravity plot shows that the drag force always assumed a minor role in determining bubble 
motion. However, in zero gravity near separation, because of the absence of buoyancy 
and the reduction in the absolute value of the dynamic force , the drag force is compara- 
ble to the other forces and is important in determining the resultant motion of the bubble. 


Effect of Surface Tension on Bubble Dynamics 

The statistical analysis (fig. 4) shows that, for the ethanol- water solution, similar 
bubble maximum-radii and lifetime characteristics were obtained in normal and zero 
gravities at comparable subcoolings. At high subcooling, this result is readily explained 
by the small role of buoyancy in normal gravity and the similar force histories, as shown 
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in figure 7. However, investigation of the force histories for low subcooling (fig. 8) re- 
veals a great difference between the normal and zero-gravity plots. Because of this 
basic change in the bubble dynamics , an alteration in the growth characteristics (average 
maximum radius or lifetime) of the bubbles could be expected. However, since the sta- 
tistical data showed that these characteristics were not affected by gravity, it is apparent 
that a change or changes took place in zero gravity in addition to the effective absence of 
gravity. 

Comparison of the ethanol-water force histories with those for water (a typical his- 
tory is shown in fig. 9) indicates that for a bubble in the ethanol-water solution the pres- 
sure force remains relatively large over the lifetime, while for a bubble in water it 
decreases rapidly from a large value near inception. The cause for this phenomenon is 
related to the shape of the bubbles. Siegel, in reference 10, suggests that for a bubble 
to be spherical the absolute magnitude of the ratio of the pressure force to the surface- 
tension force is 1. Therefore, a measure of the distortion of a bubble can be obtained 
from the deviation of this ratio from 1. As shown in figure 10, where the distortion ratio 
is plotted as a function of dimensionless time (the ratio of real time to lifetime), the 
ethanol-water bubbles were, in general, more spherical than water bubbles. Further, 
at low subcooling the ethanol-water bubbles in zero gravity were more spherical than 
bubbles generated in the same liquid in normal gravity. This condition resulted in the 
pressure force dominating bubble dynamics for this liquid at low subcooling in zero grav- 
ity and at high subcooling for both gravity levels. Since the maximum distortion of a bub- 
ble occurs during collapse, the ethanol-water bubbles may be more spherical because of 
their characteristic proportionately longer growth times (fig. 10). 

In summary, there was little change in the average maximum radius and lifetime of 
bubbles in a boiling ethanol-water solution when the gravity level was essentially re- 
moved. This condition is attributed, at high subcooling, to the dominance of the pres- 
sure force in both normal and zero gravity and, at low subcooling, to the increased im- 
portance of the pressure force in zero gravity as compared with normal gravity. The 
predominance of the pressure force is caused by the spherical shape of the bubbles. 


SUMMARY OF RESULTS 

An experimental and analytical study of the effect of an increase in viscosity in com- 
parison with the viscosity of water, on the dynamics of vapor bubbles generated on a 
flat surface yielded the following results for a sucrose solution: 

1. During the growth phase, the viscous effects on bubbles in the direction perpen- 
dicular to the heater surface are negligible. 
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2. The drag force is of consequence near separation for bubbles generated in a rela- 
tively viscous fluid in zero gravity. 

A similar study of the effect of a reduction in surface tension in comparison with that 
for water indicated the following results for an ethanol-water solution: 

3. At a particular subcooling, the average -lifetime and maximum-radii characteris- 
tics of bubbles were similar in normal and zero gravities. 

4. The dimensionless growth times were longer than for comparable test conditions 
in water. 

5. There was no necking down of the bubbles, and their contact angles remained 
relatively small at separation. 

6. The bubbles remained more spherical than those in water, resulting in a greater 
pressure force over a longer portion of the lifetime. Consequently, at high subcooling in 
both gravity levels and at low subcooling in zero gravity, the pressure force dominated 
bubble removal. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 5, 1967, 

124-09-03-01-22. 


9 


APPENDIX A 


SYMBOLS 


diameter, ft; m 

buoyancy force, lb force; N 

drag force, lb force; N 

dynamic force, lb force; N 

pressure force, lb force; N 

surface-tension force, lb force; N 

acceleration due to gravity, ft/sec ; m/sec 

gravitational constant, lb mass/(lb force)(ft/sec ); kg/(N)(m/sec ) 

constant peculiar to each liquid 

pressure, lb force/sq ft; N/sq m 

radius measured from bubble center of mass, ft; m 

radius of curvature of top surface, ft; m 

radius measured from nucleation site, ft; m 

temperature, °F; °C 

time, sec 

bubble total volume, cu ft; cu m 

bubble volume directly over base, cu ft; cu m 

velocity, ft/sec; m/sec 

velocity in r-direction, ft/sec; m/sec 

velocity in 0-direction, ft/sec; m/sec 

distance above heater surface to bubble center of mass, ft; m 
angular displacement of r, rad 

dynamic viscosity, (lb force)/(sec)(sq ft); N/(sec)(sqm) 
density, lb mass/cu ft; kg/cu m 
surface tension, lb force/ft; N/m 
dimensionless time 



cp contact angle 

Subscripts: 
b base 

l liquid 

max maximum 
S heater surface 

sat saturated conditions 

v vapor 

w water 


APPENDIX B 


SURFACE-TENSION FORMULATION 

The surface tension of water as a function of temperature was determined from 

or = 51. 8X10“ 4 (l - — Y' 2 (Bl) 

w \ 694. 4/ 

as presented in reference 11. The relation between the surface tension of water and 
some other aqueous organic solution is presented in reference 12 as 

— = 1 - K (B2) 

a w 

where K is a constant peculiar to each liquid. Inserting equation (B 1) into equation (B2) 
and solving for or results in 

cr = (1 - K) (51. 8X10 -4 ) (l — Y 2 (B3) 

\ 694.4/ 

The constant K is determined from equation (B3) by inserting a known value of a at 
some temperature. The known value of cr for the ethanol-water solution was taken di- 
rectly from reference 6, and for the sucrose solution it was obtained by extrapolating 
data from that reference. The approximations obtained from this formulation are con- 
sidered to be reasonably accurate because the liquids were far from their respective 
critical points. 
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TABLE I. - HEAT-TRANSFER DATA FOR BOILING FROM A HORIZONTAL FLAT SURFACE 



TABLE H. - COMPARISON OF MAGNITUDE OF VISCOUS AND LIQUID INERTIA TERMS 


FROM EQUATION OF MOTION (EQ. (10)) 


Angular 

Time, 

j Liquid inertia terms 

Viscous term, 

displace- 
ment of r, 

t, 

sec 

J 

PjV r (dv T / dr) 

H 

° v r 

e. 


lb mass 

kg 

lb mass 

kg 

r 2 

se 2 

deg 


(sq ft)(sec 2 ) 

(sq m)(sec ) 

(sq ft) (sec 2 ) 

p 

(sq m)(sec ) 

lb mass 

kg 







(sq ft) (sec 2 ) 

p 

(sq m)(sec ) 

10 

19X10 -5 

380 

605 

-1.23X10® 
-3. 68X10 4 
-1. 13X10 4 

-6. 00X10® 
-17.96X10 4 
-5.51X10 4 

-1.24X10® 
-3. 57X10 4 
-1.79X10 4 

-6.05X10® 
-17.42X10 4 
-8. 73X10 4 

-1. 86X10 1 
-6. 91X10" 2 
-2. 01X10" 2 

-9.07X10 1 
-33. 72xl0“ 2 
-9.81X10' 2 

30 

19x10“® 

380 

665 

-2.85X10® 
-2.51X10 4 
-6. 50x10® 

-13.91X10® 
-12.25X10 4 
-31. 72X10® 

-9. 95X10 4 
-2. 76X10 4 
-7. 49X10® 

-48.6 xlO 4 
-13.47X10 4 
-36. 55X10® 

-6.35 

-2.86X10' 2 

-1.38X10' 2 

-30.99 
-13. 96xl0“ 2 
-6. 73xl0“ 2 

60 

19x10“® 

380 

665 

-1.71X10® 
-1. 57X10 4 
-4. 87X10® 

8. 34x10® 
7. 66X10 4 
23. 76x10® 

-9. 18X10 4 
-1.23X10 4 
-5. 79X10® 

-44. 80X10 4 
-6.00X10 4 
-2.82X10® 

2.24 

2.5 XlO' 4 
-1.64X10' 2 

10.93 

12.20X10' 4 
-8.0 X10' 2 

90 

19X10“® 

380 

665 

-2. 04X10® 
-1. 08X10 4 
-9. 72X10® 

-9. 96X10® 
-5.27X10 4 
-47.4 X10® 

-1.5 X10® 
-1.46X10 4 
-6. 18X10® 

-7.3 xlO® 
7. 12X10 4 
-30. 16x10® 

9.89 

2.61X10 -2 

-1.79X10' 2 

48. 26 

12.74X10 -2 
-8. 74X10” 2 
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Figure 1. - Boiler apparatus. 
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Figure 2. -Effect of temperature on surface tension, viscosity, and spe- 
cific gravity. 





Figure 3. - Coalescence of successive bubbles in sucrose solution. 
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Figure 4. - Effect of subcooling and fluid properties on bubble maximum radius and lifetime. 
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gure 6. - Dynamics of sucrose bubble at approximately 15° F (8. 34° C) subcooling in 
normal and zero gravity. 
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